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Systematic 1/N C counting of correlators is performed to directly relate quark-gluon dynamics to 
qualitatively different hadronic states order by order. Both 2q and 4q correlators of a quanta are 
analyzed with 1/N C separation of the instanton, glueball, and in particular, two meson scattering 
states. The bare resonance pole with no mixing effects are studied with the QCD sum rules (QSR). 
The bare mass relation for large N c mesons, m p < m{^" , ^° < m^ -7-0 , is derived. The firm 
theoretical ground of the QSR at the leading 1/N C analyses is also emphasized. 
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Quantum chromodynamical (QCD) descriptions of 
hadron properties have direct relevance to understanding 
of the nonperturbative aspects of the strong interaction. 
For instance, success of the constituent quark model for 
global hadron spectra has illuminated some properties of 
chiral symmetry breaking and confinement, and provided 
the concept of constituent quarks as quasiparticles inside 
hadrons. Yet, not all hadrons are alike. In addition to 
the well-established Nambu-Goldstone bosons, there ex- 
ist some exceptional states, for example, the light scalar 
mesons (a, k, a , f ) [1], newly observed charmonia (X , 
Y, Z) 0], baryon resonance A(1405) and flavor exotic 
© + (1540) [3j. The studies of these exotic hadrons could 
provide new viewpoints beyond the simple consitituent 
quark picture, such as multi-quark components and/or 
inter- hadron-dynamics. 

The lightest scalar meson, a with / = J = 0, is a typi- 
cal example which has almost all important ingredients as 
the exotic hadrons. The a meson includes not only usual 
qq (2q) but hadronic states beyond the constituent quark 
picture: glueball, 7T7t molecule, and qqqq (4q) state with 
diquark correlation Thus it is a good laboratory to 
explore the properties and interplay of these states. Ex- 
traction of these properties has direct phenomenological 
importance to not only the hadron spectroscopy but also 
nuclear/quark- hadron mater, through the properties of 
nuclear force 0], chiral order parameter of QCD [6[, and 
diquark correlation [4(. 

The studies of the a meson, however, are not straight- 
forward. Since a can be described as admixture of several 
hadronic states, it is difficult to identify which hadronic 
states are responsible to which part of the a properties. 
Experimental information is still not satisfactory to de- 
rive the definite conclusion. Therefore, as a first step, 
it is important to theoretically clarify the properties of 
each hadronic state in the absence of mixing with other 
states 0] and, at the same time, to illuminate the role 
of interplay between these hadronic states, by examining 
the descrepancies between states with no mixing effects 
and the experimental a data with all mixing effects. 

For this purpose, we introduce the classification based 
on inverse expansion of number of color, 1/N C [1, ||, for 
the hadronic states in the correlators made of quark- 



gluon fields: U{q 2 ) = i J d 4 x e iq - x (TJ{x) Jt(0)), includ- 
ing contributions from all possible hadronic intermediate 
states with the same quantum number as J{x). One of 
the largest virtues oil/N c expansion is that it directly re- 
lates the 1/N C classifications of the quark-gluon graphs to 
the qualitative classifications of the hadronic states with 
the same quantum number, in a way that the mixing 
of these hadronic states are suppressed by higher order 
quark-gluon graphs of 1/N C . Then we can concentrate 
on the graphs for the hadronic states of our interest, sep- 
arating mixing effects from higher order of 1/N C . 

In this letter, we demonstrate this idea in the case of 
the correlators of 2q and 4q operators with the a quan- 
tum number. On the basis of the l/N c distinction, we 
can give an inductive definition of the bare "2q" and 
"4q" states, being free from the contributions of glue- 
ball, instanton, and, in particular, 7T7t scattering states, 
which are the origin of the large width ~ 500 MeV [T3| 
and large background in the a meson spectrum. A novel 
consequence of our approach is that the bare "4q" state 
can be investigated independently of 2-meson states, and 
this explicitly demonstrates the efficiency of the 1/N C 
distinction of states with the same quantum number but 
with qualitative differences. The existence and proper- 
ties of the "4q" state are dynamically studied with com- 
paring them to the "2q" state through the QCD sum 
rules (QSR) [ll|, whose theoretical ground is firm in the 
leading l/N c . We will figure out the importance of the 
"4q" component with smaller mass than "2q" case by 
150 ~ 200 MeV despite of larger number of quarks par- 
ticipating in the dynamics. This indicates the existence 
of the nontrivial correlation for the mass reduction of the 
"4q" system. 

The interpolating fields used in this work are sum- 
marized as follows: The 2q interpolating fields are de- 
scribed as = qTF^Aiq, where Dirac matrix Tm la- 
beled by M = (S,P,V,A,T) for (1^75,7^,7^75,^), 
respectively, and tf (F = 1,2,3) are the Pauli matrices 
acting ong= (u, d) T . The 4q operators with the a quan- 
tum number are given (assuming the ideal mixing for the 

a meson) by Jmm(x) = Sf=i Jm( x )Jm( x ) as products 
of meson operators (Hereafter we take the SU(2) chiral 
limit for simplicity). 
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Here we first see the 1/N C linking between quark-gluon 
dynamics and hadronic states in the case of the 2q cor- 
relators. The known facts are as follows: [!,[§]: 1) For 
quark-gluon diagrams, n-internal quark loops are sup- 
pressed by l/N™. In terms of hadron, n-meson or mul- 
tiquark with (qq) n production from Jm are suppressed 
by l/N™. 2) The disconnected diagrams with two gluon 
emission are suppressed by 1/N C . In terms of hadron, 
gg-glueball mixing is suppressed by 1/N C . 3) Instanton 
effects are suppressed by ~ e c . 4) For the meson 
properties the n-meson couplings are given as g n M — 
0(Nc~ n ^ 2 )- Here 1)~3) implies that leading diagrams 
with 0(N C ) can be naturally interpreted as the bare "2q" 
state since the 4q propagation diagrams/7T7r scattering, 
glueball, and instanton contributions do not appear at 
this order. 

Similar identification and subsequent separation are 
also possible in the case of the 4q correlator, which in- 
corporates the 4q participating diagrams from the begin- 
ning. Here we give inductive definition on the "4q" state 
following the \/N c based orthogonal conditions: a) "4q" 
can not appear in the leading N c 2q correlator, b) "4q" 
can appear in the 4q correlators even after the separation 
of 2-meson scattering states. These conditions insist that 
its dynamical origin is different from "2q" and 2-meson 
molecule states (The glueball and instanton are easily 
verified to be higher order in 1 /N c than those considered 
below, thus we will not discuss them in the following). 

Although this definition gives a convenient start- 
ing point to discuss the qualitative difference between 
hadronic states in the a meson, the studies of the "4q" 
component require systematic 1/N C arguments for the 4q 
correlators, beyond the leading 0(N 2 ) quark-gluon dia- 
grams including only 2 planar loops (FigfTJa), which are 
naturally interpreted as free 2-meson scattering and are 
irrelevant for the studies of the "4q" properties. Thus we 
must proceed to the next leading order of l/N c , 0(N C ) 
diagrams which could include 2-meson scattering, " 2q" , 
and "4q" states at this next leading order of 1/N C . 

The 0(N C ) quark-gluon diagrams in the 2-point func- 
tion could include the various hadronic contributions. 
The easiest way to classify them is to consider the over- 
lap strength of the operator J mm with hadronic states, 
involving all the elements of hadronic diagrams. We first 
classify the overlap strength of the 4q field with the 2 
meson states based on 1/N C , employing 3-point correla- 
tor among the 4q current Jmm and two separated me- 
son operators Jm' (FigdJ d-f). An explicit examina- 
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tion of quark-gluon graphs shows that the leading or- 
der diagrams are 0{N 2 ) for M — M' case, and 0(N C ) 

for M 7^ M'. Combining these facts with the overlap 
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strength of J' M with the 2q meson state \M ) is 0(N C ), 



the remaining part should be 

(0| J M m\M'M') = 0{N c )8 M m> + 0(1) 



(1) 



The first term represents the direct coupling to \MM), 
while the second term reflects that transition into 
\M'M') final state needs the additional interactions. 
This higher order counting is crucial for the separation 
of the 7T7T scattering states from 2-point correlators. 

On the other hand, the overlap strength with "2q" 
and "4q" states cannot be deduced from l/N c arguments 
only, and are assumed to be 

(0\Jmm\R) = 0(^ c 1/2 ), (R = "2q" or "4cf) (2) 

which will be assured later through the dynamical cal- 
culations. Similarily, the "4q" state will be identified by 
examining the quantitative difference of poles in 2q and 
4q correlators, which is found to be large enough to dis- 
tinguish the "4q" and "2q" states. The coupling of R to 
two mesons is estimated by 3-point function in the same 
way to obtain the meson couplings performed in Refs. Q 

— 1/2\ 

and is found to be 0(N C ). 

Now we can classify the hadronic states in the 2-point 
correlators (JmmJm'M'} based on 1/N C (See, FigQ] a-c): 
(i) If M = M', 0(N 2 ) quark-gluon graphs include only 
the free 2M scattering states in the region E > 2m m- 
Otherwise, the contributions from these quark-gluon di- 
agrams vanish, indicating absence of 2 meson scattering 
states, (ii) 0(N C ) graphs include the 2M or 2M' scatter- 
ing and possible resonance, "2q" and/or "4q". 

Note that the relations {1} in the case of M, M' ^ 
P nor A indicate that the 2-7T scattering states are not 
included up to 0(N C ) diagrams, and then the resonance 
peaks (if exist) below 2him are isolated and have zero 
width since the decay channel is absent. Therefore, now 
we can reduce the a spectrum in the 4q correlator into 
peak(s) plus continuum if we retain only diagrams up to 
0{N C ). 

This separate investigation of the 0(N 2 ) and 0(N C ) 
part of QCD dynamics enables to perform the step by 
step analyses for the 2 meson scattering, "2q", and "4q" 
spectra. In the application of the QSR, we perform the 
operator product expansion (OPE) for the correlators in 
deep Euclidean region (q 2 = — Q 2 ), then translate them, 
term by term of 1/N C , into the integral of the hadronic 
spectral function through the dispersion relation: 



nT4-Q z ) = da- 

' 7T S + 



lining (a) 



(n = 2,l) 



(3) 



FIG. 1: (Color online) Examples of the 0(iV c 2 ) and 0(N C ) 
quark-gluon diagrams for 2 and 3 point correlators. 



Now we emphasize the practical aspects of 1/N C ex- 
pansion in the application of the QSR. First, the higher 
dimension condensates in the OPE, whose values have 
been not well-known despite of their importance, can 
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be factorized into the products of known condensates, 
(qq), (G 2 ), and (qg a aGq). For example, ((qQ)(Qq)) = 
(qQ)(Qq) +J2(0\<lQ\M)(M\Qq\0) (qQ)(Qq), holds 
in leading 1/N C estimation since (qq) is 0(N C ), while 
(0\qTq\M) is 0(nI' 2 ) [|. Keeping this merit, we will de- 
duce the final 0(N C ) results from the off diagonal correla- 
tor (JvvJss)' whose leading order is 0(N C ) thus without 
the factorization violations at the 0(N C ) OPE. 

Secondly, the lowest resonance in the reduced 0{N C ) 
spectra Imll^ (s) for the "2q" and "4q" states, can be 
described as the sharp peak because of the absence of 
the decay channel. Applying the usual quark-hadron 
duality approximations to the higher excited states, 

-rh 



Vector (S=1, L=0), O(Nc) only 



Scalar (S=1, L=1), O(Nc) only 



Trlmn^(s) = \ 2 S(s - m 2 h ) + 8{s - s th )7:ImU ^(s), and 
after the Borel transformation for Eq.([3]), we can express 
the effective mass as 



m 2 (M 2 ;s th ) 



f* th ds e- s ' M \ Imn°P e (s) 
\ th ds e- s / M2 Imn°P e (s) 



(4) 



Sth can be uniquely fixed to satisfy the least sensitiv- 
ity [l6[ of the expression ((4J against the variation of M, 
since the physical peak should not depend on the artifi- 
cial expansion parameter M. This criterion is justified 
only when the peak is very narrow, and our 1/N C re- 
duction of spectra is essential for its application to allow 
the QSR framework to determine all physical parameters 
(mh, A, Sth) in self-contained way. 

In practical application of the QSR, it is essen- 
tial to reduce errors of finite order trunction of OPE 
and of the quark-hadron duality approximation. Thus 
Eq. (|4]) must be evaluated in the appropriate M 2 window 
(M^ in , M^ iax (st/ l )) to achieve the conditions: good OPE 
convergence for M min (highest dimension terms < 10% 
of whole OPE) and sufficient ground state saturation for 
M max (pole contribution > 50% of the total) [II [II [lj] . 
Without the M 2 constraint, we are often stuck with the 
pseudo-peak artifacts [l4| often seen in multiquark SRs. 
Thus we carry out OPE up to dimension 12 15] to include 
the sufficient low energy contributions which is essential 
to find the reasonable M 2 window [H, [3] . 

We summarize the numerical values used in the analy- 

— 1/2 

ses. The gauge coupling constant behaves like 0(N C ), 
and the condensates, (O) = ((qq), (a s G 2 ), (qg s crGq)) 
are 0(N C ). Here we additionaly put simple N c scal- 
ing assumptions: a s |jv c = 3a s /N c , (0)\n c = (0)N c /3. 
We use the following values with errors for the N c — 
3 case, a s (lGeV) = 0.4, (a s G 2 /ir) = (0.33 GeV) 4 , 
(qq) = -(0.25 ± 0.03 GeV) 3 , and (qg s aGq) / (qq) = 
(0.8 ± 0.1) GeV 2 . The results shown below will be ob- 
tained with the central values. We finally show (qq) and 
Wq dependence of masses in a wide range since most of 
errors come from these values. 

We start the Borel analyses from the case of the large 
N c 2q correlators for the vector meson as a reference 
and the scalar meson as the "2q" state in the a meson. 
The OPE results (up to dimension 6) are nothing but re- 
sults in the literature with employing the factorization. 
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FIG. 2: (Color online) The 0(N C ) effective mass plots for vec- 
tor and scalar mesons in the cases of various E t h values. The 
downward (upward) arrow represent the MV„ {M^ ax (sth))- 
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FIG. 3: (Color online) The effective mass plots for (JssJ vv ), 
including only 0(N C ) OPE diagrams. The large N c "2q" 
scalar meson mass is also indicated as a reference. 

Shown in Fig|2] are the effective mass for the large N c 
vector and scalar meson as functions of M 2 for various 
E t h- The downarrow (upperarrows) indicates the value 
ofM mm ( M max( s th))- Following the E th (= ^h) fixing 
criterion, we obtain Eth to 1.0 (1.4) GeV for the vector 
(scalar) meson, and determine the mass as 0.65 (1.10) 
GeV. The mass splitting ~0.45 GeV between the vec- 
tor and scalar mesons roughly coincides with the angular 
excitation energy expected from the naive consitituent 
quark picture. The reasons for slightly small value of the 
large N c p meson mass could be due to the absence of 
the factorization violations [l7| if the large N c scaling 
of condensates holds. Since our original interests are in 
the qualitative difference in large N c mesons rather than 
their absolute values, thus we will not discuss the details 
of the absolute values further. 

Now we turn to the 4q correlator results of our main 
interest. We have investigated 0(N 2 ) and 0(N C ) part of 



(JssJ ss ), (JvvJvv), and 0(N C ) part of (JvvJ S s)- We 
have checked that the typical effective masses for 0(N 2 ) 

part of (JssJss) ((■JvvJvv)) are weu above the twice of 
the large N c meson mass 2.2 (1.3) GeV, indicating that 
there is no prominent structure below the free 2-meson 
threshold, as expected from l/N c arguments. 

The "4q" state can appear from 0(N C ) part. Shown 
in FigJ3]are the effective masses deduced from (Jvv J : 



t \ 
ss / 

for E th =1.0, 1.2, and 1.4 GeV. We take the E th =1.2 
GeV case and evaluate its mass as ^0.90 GeV, which 
is obviously lower than that of the "2q" scalar meson 
case, ~1.10 GeV in large N c limit, and thus is consid- 
ered as the mass of the "4q" state. The threshold value 
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FIG. 4: (Color online) The condensate value dependence of 
masses of the large iV c mesons, "4q" (mo=0.7, 0.8, 0.9 GeV 2 ), 
"2q" scalar and vector states. 

1.2 GeV for E t h, much below the 2 scalar (vector) meson 
threshold, 2.2 (1.3) GeV, may be due to the "2q" scalar 
meson contribution, since our 0(N C ) 4q correlators can 
also include the " 2q" contribution. We have also investi- 
gated the 0(N C ) part of (JssJgs) ((JvvJyv))' an< ^ ob- 
tained the almost same mass 0.80 (0.90) GeV although 
they could suffer from the factorization violation coming 
from 0(N 2 ) OPE. These 3-independent 0(N C ) correlator 
analyses consistently suggest existence of the "4q" state 
lighter than the " 2q" state. 

Finally, we derive a conservative conclusion which does 
not depend on the details of our numerical parameters, 
especially (qq) and m\ (The dependence on the other 
parameters is relatively small). Shown in FigfJ]are "2q" 
vector, scalar meson masses, and of the "4q" mass (de- 
duced from (JvvJss)) as functions of the (qq) and m\. 
We found that the inequality m p < m\= J=0 < m^ 
holds irrespective of details of the condensate values. 

The results obtained here suggest that the a meson has 
the "4q" component, which is not generated from 7T7T in- 
teractions but from those at the quark-gluon level. Here 
we comment on Pelaez's elaborated work for the a me- 
son using the unitarized chiral perturbation with l/N c 
expansion [l8[ . He showed that the a state as the pole in 
the 7T7T scattering amplitudes of T-matrix disappears in 
large N c limit in contrast to the case of ordinary mesons 
such as the p meson. This is not contradicted with our 
results since the 7T7r-"4q" mixing is suppressed by 1/N C 



and "4q" state is not accessible from the 7T7t initial states. 
This is in sharp contrast to the 4q correlator approach 
including the "4q" state directly generated from 4q cur- 
rent. 

We conjecture that the "4q" component may play im- 
portant role as a building block of the a meson. To de- 
velop this possiblity, we plan to study the 3-point cor- 
relator for the "4q"-7T7r coupling strength. The coupling 
strength should be large since the a in nature is the broad 
resonance. If this is indeed the case, the a meson could 
be described as the 4q core clothed by the itit clouds. 
The relative importance of the " 4q" state in the a meson 
can be investigated through the coupled channel analyses 
using the effective Lagrangian including not only the it 
field but also an elementary "4q" field, whose effects are 
considered to be hidden in the parameters or regulation 
constants in the usual chiral perturbation approaches. 
This is somehow related to the recent arguments for the 
iV*(1535) resonance 0|- 

l/N c arguments developed in this work are expected 
to have the wide applications. First, the spectrum re- 
duction is applicable to the analyses of other tetraquark 
candidates. We have already obtained results that the 
effective mass in / = 2, J — channel does not show 
any stability against M 2 variation, which indicates ab- 
sence of "4q" states consistently with experiments. Sec- 
ond, the QCD sum rules which are firm in the large N c 
limit could provide useful information to models based on 
the Gauge/Gravity duality (large N c QCD <-> SUGRA), 
through the properties of the lar ge N c mesons. All these 
issues will be reported in future [20j. 
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